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Abstract 

The Fermi-Pasta-Ulam problem was one of the first computational experiments. It has stirred 
the physics community since, and resisted a simple solution for half a century. The combination 
of straightforward simulations, efficient computational schemes for finding periodic orbits, and 
analytical estimates allows us to achieve significant progress. Recent results on (/-breathers, which 
are time-periodic solutions that are localized in the space of normal modes of a lattice and maximize 
the energy at a certain mode number, are discussed, together with their relation to the Fermi- 
Pasta-Ulam problem. The localization properties of a (/-breather are characterized by intensive 
parameters, that is, energy densities and wave numbers. By using scaling arguments, (/-breather 
solutions are constructed in systems of arbitrarily large size. Frequency resonances in certain 
regions of wave number space lead to the complete derealization of (/-breathers. The relation of 
these features to the Fermi-Pasta-Ulam problem are discussed. 
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I. INTRODUCTION 



The characterization of dynamical excitations of a system is one of the fundamental 
tasks in condensed matter physics. For quantum systems the problem is obtaining the 
correct description of the many-particle wave functions. For many-particle nonlinear classical 
systems the main difficulty is relating the excitations to families of trajectories in phase space. 

Because a generic classical Hamiltonian system is non-integrable, it will evolve chaotically 
on large observation times for almost any initial condition. Hamiltonian chaos describes 
the decay of excitations on large time scales. The excitations are typically characterized 
by nonchaotic trajectories. Kolmogorov, Arnold, and Moser (KAM) showed that regular 
dynamics on invariant A-dimensional tori persists in many non-integrable systems with a 
finite number of degrees of freedom N. 1 Especially important are the simplest realizations of 
low-dimensional invariant structures - periodic orbits. Periodic orbits exist even in strongly 
chaotic systems and can densely fill the chaotic phase space volume despite having zero 
measure (just like the rational numbers fill the space of real numbers). Periodic orbits are 
the most natural class of regular trajectories which can describe excitations. An example 
of such excitations are discrete breathers. They are periodic solutions of nonlinear lattice 
systems localized in the direct space that persist in infinite lattices.-^ 

Periodic orbits can be computed numerically with high accuracy for large systems. There- 
fore computational approaches can often help in characterizing dynamical excitations, es- 
pecially when combined with analytical arguments. In the following we will show how such 
methods allow us to successfully study one of the most famous problems in computational 
and statistical physics - the Fermi- Past a- Ulam problem. 

II. THE FERMI- PASTA-UL AM PROBLEM 

In 1953 Enrico Fermi, John Pasta, and Stanislav Ulam studied the problem of energy 
redistribution between the linear normal modes of an anharmonic atomic chain and published 
the results as an internal report of the Los Alamos National Laboratory- They used one of 
the first available computers with modern architecture, and performed what can be called 
one of the first computational experiments. 

The equations of motion describe N particles with a nonlinear interaction between nearest 
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FIG. 1: The evolution of the energy of modes q = 1, 2, 3, 4 for the FPU trajectory (oscillating 
curves) and the exact periodic orbit of the g-breather (horizontal lines). The parameters are 
a = 0.25, total energy E = 0.077, and N = 32M 



neighbors, with either quadratic (the a-FPU model) 
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or cubic terms (the /5-FPU model) 
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The coordinate x n (t) corresponds to the displacement of the nth particle from its equilibrium 
position. The particle mass and harmonic spring constant are assumed to be one. 
Fixed boundary conditions Xo = xn+i = are used. The canonical transformation 

— — N 



x n (t) 



N + 



sm 



q=l 



irqn 
N+1 



(3) 



diagonalizes and solves the linear problem a = (3 = using the normal mode coordinates 
Q q (t). The mode number q — 1, . . . , N relates each of these modes to their corresponding 
normal mode frequency uo q = 2 sin (irq/2(N + 1)). The equations of motion in normal mode 
space become 



Qq + UqQq 
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for the a-FPU chain in Eq. (CQ) and 
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FIG. 2: (a) Distributions of the mode energy densities for the FPU trajectory with go = 1, N = 31, 
a = 0.33, and E = 0.32. Circles: t = 10 4 , squares: t = 10 5 , diamonds: t = 10 6 . The dashed line 
is the g-breather from (b) for comparison, (b) Distributions of the mode energy densities for the 
(/-breather with the same parameters as in (a)* 1 ^ 

for the /3-FPU chain in Eq. ([2]). The interaction coefficients 

Bqj^m = ^^(^±i±m,0 — $ q±l±m,2{N +1) ) > (6) 
± 

Cq,l,m,n = {5q±l±m±n,0 ~ ^g±Z±m±n,2(jV+l) — $q±l±rn±n-2{N+l)) (7) 
± 

give the coupling between the modes. The coupling is nonlinear, selective, and does not 
have a characteristic interaction range (or distance) in normal mode space. Thus, even 
very distant modes in normal mode space are coupled. In the absence of this interaction 
there exist iV integrals of motion which correspond to the energies of the normal modes 
E q = (Ql + ^ q Ql)/2. 

Fermi, Pasta, and Ulam considered the evolution of the nonlinear chain with the single 
mode q = 1 with the lowest frequency initially excited.- They used standard numerical 
schemes to approximate the differential equations and monitored the state of the system 
for times which are orders of magnitude larger than the largest periods of oscillation of the 
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normal modes. It was expected that after some time other modes would become excited 
as well, and the energy would eventually be equally distributed over the entire spectrum. 
That was the way they anticipated a transition to thermal equilibrium. To their big surprise 
they observed the opposite result (see Fig. [[]). The energy stayed localized in several low 
frequency modes for times that are orders of magnitude larger than the typical oscillation 
periods of these normal modes. Recurrences of almost all the energy into the initially excited 
mode were observed as well. Later computational studies revealed that there exist energy 
thresholds and system size thresholds above which equipartition was observed on relatively 
short time scales.-^ Many questions arise. Why is the energy localized in a few modes below 
the thresholds? Why does the energy spread quickly into the entire spectrum above the 
thresholds? How do the threshold values scale with the parameters of the system? 

Figure [2] shows the mode energy density distribution e q = E q /N for an initial energy 
that is larger than the one considered in Ref. 5;, but yet well below the above mentioned 
thresholds. The distributions are plotted on a logarithmic scale and for times t = 10 , 10 5 , 
and 10 6 . Because the mode energies are time-dependent, the densities were averaged over 
a time interval of At = 10 4 . We observe exponential localization of the energy distribution 
in normal mode space, which was first noted by Galgani and Scotti.- Such a distribution is 
formed on a very short time scale 7~i,and possesses a core — the few modes which are strongly 
excited — and a tail of exponentially weakly excited modes. Also shown in Fig. [2] is a slow 
resonant pumping of energy from the core of the distribution into its tails for later times. 
This process ultimately brings the system to equipartition, but on the much larger time 
scale T2 ^> T\. Fermi, Pasta, and Ularn^ performed their numerical studies at lower energies, 
and therefore r 2 was even larger in their case. If the initial energy is further increased, a 
threshold value may be reached for which the time scale of the core-tail pumping becomes 
comparable with the transition time to the exponential distribution r 2 pa ti. 6,7 Below these 
threshold values the time scales differ by orders of magnitudes, and the evolution on times 
ri < f < T2 corresponds to the observations of Ref. la. Fermi, Pasta, and Ulam viewed 
their observations as a problem of (perhaps complete) the absence of equipartition (that is, 
ergodicity) - a problem now known as the FPU problem. From a contemporary perspective, 
the FPU problem can be reformulated as the following set of questions: 

• Why are there at least two qualitatively different time scales t\ and r 2 ? 
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• Are there low- dimensional invariant manifolds in phase space which approximate the 
dynamics of the FPU trajectories for times less than T2? 

• Can we estimate the localization length of the exponentially localized energy distribu- 
tion in normal mode space? 

• Is the localization length related to the time scales? 

• Does the FPU observation hold in the limit of infinite systems, as well as in two- and 
three-dimensional systems? 

There is a long history of studies of the FPU problem,-^^ which, even if not directly 
answering all of these questions, defined new fields of research in nonlinear dynamics, chaos 
theory, and soliton theory. 

III. PERIODIC ORBITS q-BREATHERS 

From the results of Sec. [TT] we may conjecture that the exponentially localized state 
observed on the FPU trajectory is a long-lived excitation, which is destroyed after time t?,- 
Therefore (see discussion in Sec. [I]) we anticipate that this excitation is close to an exact 
periodic orbit for times i<T 2 . Then the observed dynamics will be almost regular. Such a 
periodic orbit will correspond to a slightly deformed normal mode. There exists a rigorous 
method of constructing such a periodic orbit, starting from the linear limit a = (3 = 0.— ^ In 
this limit the periodic orbit is given by exciting a single mode with mode number q with an 
energy E. With these conditions we isolate a unique periodic orbit in the multi-dimensional 
phase space of the system. 

The linear spectrum of the FPU chain obeys the non-resonance condition^ 

uj q ^ nu qo (8) 

for all integer n and q ^ go- This condition allows us to apply a theorem due to Lyapunov,— 
which states that the original periodic orbit can be continued into the domain of nonzero 
nonlinearity for fixed energy— Continuation means that close to the periodic orbit of the 
linear system we will find a new periodic orbit in the phase space of the nonlinear system 
with identical energy E. This new periodic orbit, being a closed loop in phase space, can 
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be viewed as a continuous deformation of the periodic orbit of the linear system, as the 
nonlinearity is increased. 

The numerical methods of constructing g-breathers^ 1 ^ use the fact that a periodic orbit 
corresponds to a fixed point in the generalized Poincare maps (sections) in phase space. Let 
us explain one of the possible realizations of such a numerical scheme in more detail. 

We consider initial conditions for which all normal mode velocities Q q (t = 0) = 0. Then 
the kinetic energy at t = is exactly zero. Note that the kinetic energy is given by the 
one-half of the sum over all squared mode velocities, and the total energy is conserved 
during the integration. We construct a iV-dimensional vector where each component is the 
amplitude of a given mode. This vector represents the configuration space of the system. 
Then we integrate the equations of motion numerically, until the velocity Q qo crosses zero 
two times (because an oscillator needs two turning points to return to its initial state). At 
that moment we again measure all mode amplitudes, and obtain a second (final) vector. If 
the initial and final vectors coincide, the potential energies coincide as well. Therefore the 
kinetic energy at the final time will be zero, and our initial condition is located on a periodic 
orbit. The period of the orbit can be determined by computing the time we needed to reach 
the final vector. 

What if the final vector does not coincide with the initial one? Then we compute the 
difference vector between the final and initial vectors. It will also have N components. 
Our task is to slightly vary the initial vector components such that all the components of 
the difference vector shrink to zero. In other words, we are trying to zero N functions 
(the components of the difference vector) by varying iV variables (the components of the 
initial vector). Generalized Newton- Raphson methods will do the job if the initial vector 
is suitably close to a desired solution. Because the final vector is the outcome of a map in 
configurational space, we can also say that we are looking for a fixed point of that map. If 
we want to find a fixed point for a given total energy (or some other constraint) we have to 
limit the variations of the initial vector accordingly. For practical reasons it is convenient 
to first fix the amplitude of mode Q qo to a given value a to find the corresponding periodic 
orbit, calculate its energy, and then to vary a such that the desired energy value is realized. 

Once a periodic orbit is obtained, we also compute its stability properties. A periodic 
orbit is said to be (marginally) stable if a small perturbation of the orbit in phase space does 
not grow in time; otherwise it is unstable. The initial smallness of the perturbation allows 
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us to perform a linearization of the phase space flow around the periodic orbit (see Ref. Q 
for details). In practice it amounts to the following procedure. Take a reference point on 
the periodic orbit. Consider 2N small perturbations in the 2N different directions in phase 
space. Integrate each of these perturbations over one period of the orbit, and obtain a new 
(still quite small) perturbation. The components of these final perturbations define a Floquet 
matrix, which maps any initial perturbation into a final one. The matrix is symplectic.— In 
a nutshell, if A is an eigenvalue of a symplectic matrix, so are 1/A, A* and 1/A*. Compute the 
eigenvalues of that matrix. Any eigenvalue which is located on the unit circle in the complex 
plane corresponds to a perturbation that does not grow in time. Eigenvalues outside the 
unit circle correspond to perturbations which grow exponentially fast in time. Note that 
each eigenvalue outside the unit circle has a symmetry-related partner inside the unit circle 
due to the symplectic properties of the Floquet matrix. The interested reader may consult 



Ref. 



171 . where details of the ways to compute and analyze periodic orbits are discussed. 



A. Results for the a-FPU model, low frequencies 

The dynamics of the continued periodic orbit incorporates many normal modes, because 
nonlinearity induces mode-mode interactions. The numerical procedure described in Sec. ITU 
was used and new periodic orbits were found. The amplitude and energy distributions of the 
normal modes in such a periodic orbit turn out to be exponentially localized (see Figs. [Hand 
[2](b) ) . We call these solutions g-breathers; they are periodic in time and localized in normal 
mode space. Note that the shape of the energy distribution of g-breathers is remarkably 
close to the one of the FPU-trajectory for t < r 2 and identical parameters (energy E and 
seed mode number g ). 

Any numerical evaluation should, if possible, be accompanied by analytical considera- 
tions. For g-breathers that is possible with the help of standard perturbation approaches. 
For the case of low-frequency modes the energy distribution in the g-breather has been 
estimated to be^ 



E =<y 2n - 2 n 2 E = "V^gol" Z ±J (Q) 



n A /E^(A + l) 3 / 2 
7T 2 q 2 

The g-breathers are localized in normal mode space if mode energies shrink with increasing 
distance from go. From Eq. (J9j) it follows that g-breathers are exponentially localized if 7 < 1: 
E q ~ e~ q ^, where £ is the localization length. For 7 = 1 Eq. ([9]) predicts the parameter 
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FIG. 3: Stable g-breathers for a = 0.025 and a = 0.25, E = 0.077, N = 32, and q = 1. The 
dashed line is the result of Eq. ([9]) .— 

values for which a g-breather delocalizes. If this condition is reached upon variation of a 
parameter, we say that a derealization threshold has been reached. If we vary a second 
parameter, we can keep 7 unchanged. In this way we can vary all parameters, keep 7 
constant, and obtain scaling relations between all the parameters. The inverse localization 
length £ _1 is proportional to q^ 1 In 7. When 7 = 1 the g-breather delocalizes, which implies 
T2 ~ T\. These are, among many others, predictions which can be tested by comparing to q- 
breather solutions and FPU-trajectories, which are obtained using the previously described 
computational methods. 

The numerical results and the estimate in Eq. are shown in Fig. [3J The g-breather 
solutions are linearly stable.— 1 ^ The agreement between the computations and theory is 
very good, especially for strongly localized energy profiles. 

Note that in real space the energy distributions are delocalized both for the FPU and for 
the g-breather trajectories, and do not show remarkable features.— 1 ^ 
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FIG. 4: (color online) The profiles of the g-breather solutions for N = 359, (3 = 1, E = 0.346, 
and qo = 4, 8, and 24. The symbols denote the mode energies at the time when the coordinates 
x n (t) = 0. The vertical bars denote the range of mode energy values realized during one period 
of the g-breather. The dashed lines are the result if Eq. (|10p with go = 24, 8, and 4 from top to 
bottom.— 

B. Results for /3-FPU, low frequencies 

If we apply the perturbation methods we have mentioned to the /3-FPU problem, we 
arrive at the result for the energy distribution in the g-breather: 14 

&q=(2n+l)q — A Jii qo , A - — ~ 2 ■ (,-W) 

07T g 

The necessary condition for the localization of the g-breather is A < 1. The inverse localiza- 
tion length £ _1 is proportional to (2qo)^ 1 In A. If A = 1, the g-breather delocalizes, and thus 
we expect that ti « t\. The numerical and analytical results ( fTUI) are presented in Fig. HI 
For q = 4 the asymptotic result for A fails, although the exponential profile of the mode 
energies still holds. 

The g-breathers become unstable if 6/3E(N + 1)/tt 2 > 1,— but as long as they remain 
strongly localized, the FPU trajectories do not delocalize quickly— Remarkably the FPU 
trajectory shows a change from regular to weakly chaotic dynamics in its core when increas- 
ing the initial energy E above a certain value.— This value corresponds to the condition for 
changing the stability of the g-breather. It is well known that for nonintegrable Hamiltonian 
systems, chaotic trajectories appear close to an unstable periodic orbit. 19 
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FIG. 5: Energy density distribution for the FPU trajectory (squares, t = 10 5 ) and g-breather 
(circles) for N = 47, E = 4.7, and q = 47. (a) a = 0.25 and (b) (3 = 0.25. 12 

Why does the instability of the g-breather, and the chaotic dynamics close to it, not 
destroy the exponential localization of the mode energy profile? The answer might be that 
a g-breather which is well localized in normal mode space does not allow for a resonant 
interaction between core and tail modes, because if such a resonance was present, it would 
lead to a derealization of the periodic orbit itself. Instabilities in the core of the distribution, 
that is, resonances between modes inside the core, are what is left. They can lead to a chaotic 
dynamics inside the core, but not to a derealization, due to the absence of resonances 
between the core and the tail modes. 

C. Results for high frequencies 

Localization is observed also for high frequency seed modes. In Fig. [5] the energy distri- 
butions for the FPU and g-breather trajectories are shown.— The only difference from low 
frequencies is that for the a-FPU model high and low frequency modes become excited in 
pairs, in agreement with the predictions of perturbation theory.— 
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FIG. 6: Mode energy distributions for a g-breather in a two-dimensional FPU lattice with N = 
32 x 32, E = 1.5, q = (3, 3), and (3 = 0.53 

IV. GENERALIZATION TO TWO AND THREE DIMENSIONS 

The condition for applying Lyapunov's theorem to the continuation of a g-breather is the 
absence of resonances (see Eq. (JSJ)). A finite system has a discrete spectrum and thus Lya- 
punov's theorem is applicable. A generalization of the results in Sec. IHII has been performed 
by considering finite lattices with spatial dimensions d = 2 and 3. Instead of mode numbers 
we now have mode vectors with d integer components. The main properties of the solutions 
for g-breather solutions do not change. 2 - In Fig. [6] an example of a g-breather solution in 
d — 2 is shown. 

For the /3-FPU model we can generalize the perturbation theory and predict the energy 
distribution^ . 

TP — \2njp \ _ 3/3-Kq A^ 

^(2n+l)q - *d ^qo> A d ~ 2 2+d 7r 2| qo |2 " 

The numerically computed g-breather solutions show quantitative agreement with these 
estimates.— 

V. SCALING AND TRANSITION TO MACROSCOPIC SYSTEMS 

The perturbation theory results in Eqs. fl9|)- f|TT|) show that after replacing the extensive 
parameters (the energy E and the mode number q) by intensive ones (the energy density 
e = E/N and wave number k = nq/(N + 1)), the localization length becomes independent 
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FIG. 7: (a) Dependence of A on the size of the chain for the g-breather with energy density 
e = 4 x 10~ 4 (lower curve) and 2 x 10 -3 (upper curve) and (3 = 1. The open circles correspond 
to N + 1 = r(No + 1) with N$ = 15, where r is an integer, (b) The dependence of A on the wave 
number ko for the corresponding data from (a)*2i 

of the system size N. Consequently, we expect that g-breathers will persist in the limit of an 
infinitely large macroscopic system. Consider a g-breather solution Q q (t) for a finite chain 
of size N and another chain of size N + 1 = r(N + 1), where r = 2, 3, 4, ... It follows that 



Qg(t) 



VrQ q (t), q = rq 



0. 



(12) 



q ^ rq 

is a solution for the larger chain.— 

By increasing r to infinity we will obtain solutions for infinitely large macroscopic systems. 
The scaling procedure is easily generalized to d = 2 and 3 as well as to free and periodic 
boundary conditions.— 

Suppose that g-breathers exist in an infinitely large chain. Then Eq. (jTUl) with k and e 
instead of q and E gives 



In* 



Sk = (j^- - l) lnVA + ln£ fco , 



3/3 gfc 
II' 



(13) 

We see that the analytical result for the mode energy profile is independent of system size 
if intensive parameters are used. 

Consider a finite chain with initial length N = 15 and seed wave number k = it/ 15. We 
increase the system size and numerically compute g-breather solutions for the seed mode 
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number go which is the closest to ko(N + 1)/tt. The solution is then used to approximate 
the coefficient A by the ratio e^ qo /e^ qo (see Fig. [3(a)). The results confirm the independence 
of the approximate value of A if N + 1 = r(N + 1). For other values of N we probe A with 
different seed wave numbers ko in the vicinity of k . The continuous behavior of A(fco) (see 
Fig.fT^b)), irrespective of system size N, confirms these conclusions for macroscopic systems. 

If the energy density e is fixed to a given value, it follows from the geometric series in 
Eq. fflOl) that Sk = (1 — A)e. We use Eq. f JT3|) and calculate the slope S of the energy density 
distribution (on a logarithmic scale):— 

a I n- rr ^1 + 4z/ 4 /^ - 1 3/5 

5 = ^ lnv ^ ^ = 2^Ao 2 ' " T e - (14) 

The absolute value of S equals the inverse localization length £ _1 = |,f>|. The slope S 
depends on the seed wave number ko and on the effective nonlinearity parameter u, which 
is determined by the product of the energy density e and the nonlinearity strength j3. The 
localization length diverges if A;o — ^ and takes a minimum value of £ m i n ~ Z//0.7432 at 
ko = fcmin ~ 2.577^. 21 The localization length for ko = k m \ n also increases for increasing 
v. In the limit k ^> v we estimate S ~ — 2/fc ln(fc /^), and for k v we obtain 
S ~ —ko/ (2i/ 2 ).— It follows that for a given energy density there exists a seed wave number 
k m in for which the g-breather is most strongly localized. Moreover, g-breathers tend to 
delocalize for k -C k min , that is, in the limit of long wavelength. These analytical predictions 
may now be compared to careful numerical computations of g-breathers. 

The analytical and numerical results are compared in Fig. [HI The numerical results 
confirm that the localization length does not depend on the system size. We also observe the 
minimum S(k Q ) whose depth and location vary with the energy as expected. A systematic 
mismatch between the theory and the numerical results for small ko is caused by higher 
order corrections to the perturbation theory. The g-breathers from Fig. H] for g = 4, 8 and 
24 correspond to the most left symbol on the middle curve in Fig. [HJ the minimum on that 
curve, and a point to the right of it, respectively. 

It follows from Eq. <HM that the quantity S m (z) = vS depends on a single variable, 
namely the dimensionless seed wave number z = ko/v. It implies that knowing this single 
master slope function is sufficient to predict the localization property of any g-breather, at 
any energy, any seed wave number, etc. In Fig. [9]the dependence of S m (z) is shown together 
with the rescaled numerical data from Fig. [HJ 
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FIG. 8: The slope S{k ) for /3 = 1 and e = 1.57 x 1(T 3 , 9.6 x 1(T 4 , and 6.08 x 10" 4 (dashed lines 
from top to bottom). The symbols and connecting lines are the results of the numerical solutions 
of g-breathers for N = 149 and N = 359 M 
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FIG. 9: S m {z) (dashed line) with z = ko/u. The symbols and connecting lines are the rescaled 
data from the numerical g-breather solutions in Fig. [HJ 21 

Despite the systematic mismatch between the numerical data and the theoretical predic- 
tion for small z = ho/ft, the scaled numerical data collapse to a single curve. Therefore, 
the theoretically predicted scaling law in Eq. (fT4"j) is well confirmed. Analogous results have 
been obtained for the a-FPU model, and for both FPU models for seed wave numbers with 
frequencies near the upper edge of the spectrum uj q . 
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VI. DISCUSSION 



The reason for the derealization of g-breathers at the edges of the linear spectrum are 
resonances between normal mode frequencies: u qo ~ nu nqo for low frequency modes and 
cu qo cu q for high frequency modes. Thus it is expected that the excited mode with a wave 
number close to an edge of the spectrum will decay into many other modes, and the two 
time scales r 2 ~ r±. These time scales must diverge in the limit of small frequencies, due 
to conservation of total mechanical momentum. Hence the time scale r S i m of the simulation 
is important. Suppose we excite a normal mode that is close to a well-localized g-breather 
when r 2 3> T\. If r sim < r 2 , the excitation of this normal mode will quickly spread into a 
packet, but stay localized in normal mode space for all observation times (this localization 
is the FPU problem). If r sim > r 2 , the long time of nonequipartition will be eventually 
replaced by a thermalized state. 

Suppose we excite a low frequency normal mode for which the corresponding g-breather 
is delocalized, and r 2 ~ T\. If r sim < r 2 , the normal mode will not decay at all during the 
observation time. The simulation will thus recover the linear dynamics of the bare modes as 
if the nonlinearity were absent. If r sim > r 2 , the normal mode will begin to decay into other 
modes, but there will be no intermediate state in which only a few other modes become 
excited; rather the entire available mode space will be excited at once. That means that 
these low frequency modes will behave like linear modes for short times and fully chaotically 
for larger times, without any intermediate regime of exciting only a few other modes. 

The normal mode frequencies ui q are periodic in the wave number k. That period defines 
the irreducible part of the wave number space, which therefore has finite width (for the chain 
it is 7r) . Localization in wave number space is meaningful only when the localization length 
£ < 71. Thus, from our previous results we conclude that all g-breathers which are closer than 
some critical distance Ak of the seed wave number /c from the spectrum edge delocalize. 
For non-zero v 2 ~ f3e Ako 7^ 0. Thus, there exists a range of k for which the excitation of 
the normal mode decays and the normal mode picture may become ill-defined. The range 
Ak is an increasing function of u, and for some critical value v Ak as ir, which means that 
normal modes do not characterize the dynamics of the system. That corresponds to the 
regime of strong interaction between modes. For small v this strong interaction occurs near 
the edges of the spectrum only. 
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Let us discuss the results for finite systems. A given v corresponds to a finite A&q. A 
finite system implies a discrete set of wave numbers for the normal modes. The thermal 
conductivity of anharmonic acoustic lattices is governed by the dynamics of low frequency 
waves. To numerically observe a conductivity different from the a free propagating (ballistic) 
one, we have to resolve the regime of strong interactions. In particular, we have to choose a 
system size greater than N c ~ l/(Ako). The theory we have presented allows us to obtain 
relevant quantitative estimates, and the computational results show that these considerations 
are valid. 

Finally, let us discuss the resonance peaks in the tails of the distributions in Fig. [21 which 
are observed in the FPU trajectory and also in the g-breather. A quantitative explanation 
of the origin and position of these peaks was given recently— They are generated by the 
closeness to the resonance in Eq. ([8]) for certain mode numbers. It is easy to check that the 
peaks appear in the upper half of the frequency spectrum and for low-frequency seed modes 
only. To explain the appearance of these peaks we have to consider the dispersion relation 
more accurately— Numerical experiments show that these resonant peaks trigger the tran- 
sition to equipartition in the FPU trajectory— The nature of the pathway to equipartition 
for high-frequency seed modes, for which these peaks are absent, remains a puzzle. 

There are many other questions that still wait for answers. We have to explore other 
nonlinear models and test which of the results we have discussed will be generic and which 
not. We need to obtain reliable estimates on the dependence of T\$ on the parameters and 
compare with the results for g-breathers. We have to determine whether g-breather-like 
excitations are spontaneously generated at thermal equilibrium and to characterize their 
statistical properties. This list could be continued, but there must be an end to everything, 
and our story finishes here. 
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